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ABSTRACT 

This paper presents a mathematical model characterizing the behavior of a simple 
amplifier using a FeFET. The model is based on empirical data and incorporates several 
variables that affect the output, including frequency, load resistance, and gate-to-source 
voltage. Since the amplifier is the basis of many circuit configurations, a mathematical 
model that describes the behavior of a FeFET -based amplifier will help in the integration 
of FeFETs into many other circuits. 
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INTRODUCTION 

The use of ferroelectric field-effect transistors (FeFETs) to create simple amplifiers has 
not been extensively studied. Moreover, the FeFET’s unique characteristics of hysteresis 
and nonlinearity result in an amplifier that displays properties that are different from 
those of a MOSFET-based amplifier [1], The FeFET amplifier’s behavior has never been 
modeled. Thus, a mathematical-based model that describes a FeFET-based amplifier was 
created. The model is based on empirical data and incorporates several variables that 
affect the output, including frequency, load resistance, and gate-to-source voltage. The 
amplifier is the basis of many circuit configurations, and creating an easy-to-use model 
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that illustrates the behavior of the FeFET amplifier will greatly aid in the integration of 
FeFETs into many other circuits. 

FeFET AMPLIFIER CIRCUIT CONFIGURATION 
A FeFET simple amplifier is built using a FeFET and a resistor. Figure 1 illustrates the 
FeFET amplifier circuit configuration. As can be seen from the figure, the resistor, R , is 
connected to the source of the FeFET. The input voltage, V in , is provided from the gate 
of the transistor to ground. Thus, V in includes both Vgs and the voltage across the 
resistor. Vdd is the voltage from the drain of the FeFET to ground, thereby including 
both Yds and the voltage across the resistor. The output voltage, V out , is measured at the 
source of the transistor and above the resistor. 



Figure 1 FeFET amplifier circuit configuration 

DESCRIPTION OF THE MODEL 

The model was created in Microsoft Excel and is based on the equations [2] for the drain 
current, Id, that are derived from the Fermi-Dirac equation. When the gate voltage is on, 
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( 1 ) 


Id = 


Idsat — B* Log(ti) 

e k{Vp-Vgs) + | 


When the gate voltage is off, 


Id - 


Idsat - B * Log(t\) 

e k(Vss ~ Vp) + 1 


( 2 ) 


In the two previous equations, Id is the drain current, Idsat is the drain current when the 
FeFET is in saturation, B is the drain current decay coefficient, V gs is the gate-to-source 
voltage Vgs, V p is a polarization voltage, which is the gate voltage at which half the 
saturation current is reached, k is a constant defining the rate of change of the function, 
and ti is the elapsed time in seconds since the last poling. More needs to be said about 
the polarization voltage, V p . Each FeFET has two different E p ’s, one is for the 
negatively-poled drain current, while the other is for the positively-poled current. This 
provides the hysteresis characteristic of the FeFET’ s I-V relationship. 

The model allows the user to enter the desired value for several of the parameters 
that affect the output. Thus, the user can define values for Vdd, load resistance ( Ri oa d ), 
input frequency, amplitude of V in ( V Uhum ,>), offset of V in ( Vi n , 0 jjsei), phase shift, and 
frequency coefficient. V in , offset was included in order for the model to be able to account 
for offsets in the input voltage. The phase shifting parameter phase is not the precise 
angular phase shift of the signal; rather it is a tunable parameter that is used to 
figuratively represent the amount of phase shift desired. Therefore, it is not a measurable 
parameter; it is an adjustable, user-specified parameter. The frequency coefficient is used 
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to regulate the effect of frequency. Other parameters incorporated into the model include 
positive and negative V p ( V p , pos and V p , ne g, respectively), k, which is the constant given in 
Equations (1) and (2), and maximum I D ,- Iomax was needed for the calculation of Id- The 
parameters V p , pos , V p , neg , k, and lomax are determined through experimentation. 

EQUATIONS OF VARIABLES 

Besides the parameters defined by the user and those obtained by experimentation, the 
model also includes variables that are determined using simple mathematical equations. 
These parameters are V uh Vj n , phased, Polarization 1 , Polarization2, Id P os, iDneg, Id, Id, fa, and 
V out . Vi n is the value of the input voltage that takes into account both the amplitude of the 
input voltage and an offset, if given. The value of the input voltage that includes V in ,am P , 
Vin.offset, and the phase shift is V in , phased- Polarizationl and Polarization2 are parameters 
that reflect the polarization property of ferroelectric material. Id PO s is the drain current 
found using positive V p , V p , pos , when the gate voltage is on, while Ioneg is the drain current 
found using negative V p , V p , neg , when the gate voltage is on. I D is the final value of the 
drain current. Id, fa is the frequency adjusted value of Id- Thus, Id, fa adjusts Id based on 
the input frequency and the frequency coefficient. Finally, V out is the output voltage. The 
equation for each parameter is given in the following set of equations. 

V in = V i nj amp* Sin(radian Value) + Vin.offset (3) 

Vin, phased = V m ,am P *sin(radian value + phase shift) + V in , offset (4) 

At radian value 0, 

Polarizationl = 0 (5) 

At all other test points, 
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Polarization 1 

1 If V in> V p,pos (6) 

= < 

the value of Polarizationl at the preceding test point if Vi n <V p ,pos 

( 7 ) 

At radian value 0, 

Polarization2 = 0 (8) 

At all other test points, 

Polarization2 

0 if V in<V Pi n C g (9) 

= •< 

the value of Polarizationl at this test point if Vi n >V Pjne g 

( 10 ) 

iDpos = iDmax / + 1) (1 1) 

iDneg = IDmax / ~ + 1) (12) 

Id 

r 

Id P os if Polarization2 =1 (13) 

= •< 

IDneg if Polarization2 = 0 (14) 

Id ; fa = lD*frequency coefficient H Togiofrequency (15) 

Vout = I D ,FA*Rload (16) 

The Vi„ and V in , phased equations are based on the sinusoidal wave equation [3] 
given by 

x(t) = A*sin(cot + 6) + D, (17) 

where A is the amplitude of the signal, co is the radian frequency in radians per second, t 
is time in seconds, Q is the phase angle in radians, and D is the DC offset. In the equation 
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for V in , Vin, am p is the amplitude, radian value is cot, and V in , offset is the DC offset. There is 
no phase shift in V in ; however, this variable exists in the equation for V inii phased- In order to 
fully explain the term radian value in the equations for V in and V in , phased, the method of 
measurement used in the model must first be discussed. The model measures the 
parameters at specific test points. The value of these test points are the nonnegative 
numbers. The degree equivalent of each test point was chosen to be (360/30) times the 
value of the test point. This relationship is given as follows 

degree value = test point * (360/30) (18) 

The radian equivalent of the degree value of each point is utilized throughout all the 
equations for the variables. Since cot is equivalent to 2nft, where / is frequency, and t in 
this model is a multiple of the period 1/f ( t=n/f where n is a nonnegative number), then 
cot can be reduced to 2k n. However, degree value can be rewritten as 360*n/30, whose 
radian equivalent is given by 2nn/30. Thus, cot in Equation (17) is replaced with 2nn/30 
in the equations for V in and V in , phased, and the term 2nn/30 is denoted by radian value. 

The equation for each of the other variables is derived from the basic properties of 
ferroelectric material. The polarization variables, Polarizationl and Polarization2, are 
initially set to zero, then altered based on the value of V in at the current test point. If Vt„ is 
greater than V p , pos , then Polarizationl is set to 1, otherwise it is set to its value at the 
preceding test point. Polarization2 is set to 0 if V in is less than V P: „ eg . If V in is greater 
than V pMe g, Polarization2 takes the value of Polarizationl at the current test point. The 
equations for I Dpos and Ioneg are derived from the Fenni-Dirac-based equation for I D when 
the gate voltage is on, where V p in Equation (1) is replaced with V PiP0S for I Dpos and V p ,„ eg 
for iDneg, Idsat is replaced with Inmax, and V gs is replaced with V in , phased- The term 
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B*log(ti) in Equation (1) is ignored in this model with little loss of accuracy. Id equals 
In,, „s if Polarization2 is 1, but if Polarization2 is 0, then Id equals loneg- V ou t is 
determined simply using Ohm’s Law with the load resistance and the frequency-adjusted 
value of Id- This set of equations along with the user-defined and experimentation-based 
parameters provide a very good approximation of the output of a FeFET simple amplifier. 

MODELED AND MEASURED DATA ANALYSIS 
The accuracy and efficiency of the model created was examined by comparing the output 
of the model with empirical data for specific test cases. It is important to note that this 
model is the first mathematical-based computer model to characterize the behavior of a 
FeFET simple amplifier. Thus, this model is a first attempt and can be considered a 
preliminary version of more sophisticated models that can be created in the future. 

Given the user-defined and experimental-based parameters, the model outputs a 
plot of V in and V out with respect to time. As was aforementioned, time is given as a 
multiple of the period. Several representative test cases were examined. For each test 
case, two sets of parameter values are included: one set leads to a model output where 
the amplitude of V ou , V out , am p, is as close to the oscilloscope value as the model allows, 
while the other set results in a waveform for V ou t that is very similar in shape to that of the 
oscilloscope. For each set of parameter values, a table of the values and a plot of the 
modeled output are given. The modeled plots are compared to the oscilloscope plot. For 
all the test cases, the sum of Vi n , amp and V in , offset never exceeds 8V to ensure that the 
FeFET is not burned. 
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Test Case I 


For the first test case, Vdd was set to 0.709V, the amplitude of V m was 6V, the frequency 
was 100Hz, and the load resistance was selected to be 35kfl No offset was included in 
the input signal. The following table lists the values of the user-defined parameters that 
resulted in V ou t,amp that is greater than that of the oscilloscope, but as small as the model 
allowed. 

Table 1 Initial parameter values for Test Case I 


Parameter 

Value 

Vdd 

0.709V 

Rload 

35 tel 

Frequency 

100Hz 

V in, amp 

6V 

V in, offset 

OV 

Phase 

1.5 

iDmax 

6.00E-04A 

V p,pos 

2.5V 

Vp.neg 

-0.5V 

k 

0.75 

Frequency Coefficient 

0.007 


The parameter phase was set to 1.5 in order to reproduce the phase shift seen in 
the oscilloscope output of the empirical measurement, which is shown in Figure 2. It was 
observed that a positive value for phase caused the output signal to be shifted to the left, 
while a negative value for phase led to a rightward shift of the output signal. Thus, in 
this test case, the output signal leads the input signal. The value of the parameter k was 
selected after experimenting and noting the effect of changing that parameter. For this 
test case, decreasing k below 1 removed the clipping effect at the peaks of the output 
signal, and further decreasing k below 0.8 removed the clipping effect at the troughs of 
the output signal. A very small value for frequency coefficient was chosen to ensure that 




V out, amp was as small as possible. Decreasing frequency coefficient decreases Id, fa, as 
seen in Equation (15), which leads in turn to a decrease in V ou t ■ The values of /b„, 
V P:POS , and Vp.neg were determined from previous experiments and will remain unchanged 
throughout all the test cases. Figures 2 and 3 show the oscilloscope output and the 
modeled output of the signals, respectively. 


Vi, 


Vout 



Figure 2 Oscilloscope output for Test Case I 


FeFET Amplifier Model 
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Figure 3 Modeled output for Test Case I based on initial parameter values 


In this and all upcoming oscilloscope plots, it must be noted that the input and 
output waveforms are not plotted on the same coordinate system, despite their sharing the 
same axes. For example, in Figure 2, the coordinate system of V !n is partitioned into 
divisions that are 2V each, whereas V ou i s coordinate system is divided into 5mV per 
division blocks. Therefore, Vi n , amp is 6V, while V out ,amp is only 1 1.8mV. 

The maximum value of V ou , as measured by the oscilloscope is 1 1.8mV, which is 
significantly larger than the amplitude of the model’s V ou t, as seen in Figure 3. However, 
this is the smallest value of V ou t,amp that still displays the shape of the signal. Comparing 
Figures 2 and 3, one might not believe that the model accurately illustrates the behavior 
of the FeFET amplifier, but increasing frequency coefficient to 0.15 and decreasing k to 
0.3 produces the plot shown in Figure 4, which closely resembles the oscilloscope output. 
The given value of the frequency coefficient was found to produce the desired amount of 
frequency effect. At 0.3, k was small enough to model the jumps in the peaks and 
troughs visible in the oscilloscope output of this test case. Table 2 lists the new 
parameter values. 
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Figure 4 Modeled output for Test Case I based on new parameter values 
Table 2 New parameter values for Test Case I 


Parameter 

Value 

Vdd 

0.709V 

Rload 

35k^ 

Frequency 

100Hz 

V in, amp 

6V 

V in, offset 

OV 

Phase 

1.5 

iDmax 

6.00E-04A 

V p,pos 

2.5V 

Vp, ne g 

-0.5V 

k 

0.3 

Frequency Coefficient 

0.15 


Comparing the modeled waveforms of Figure 4 with the measured waveforms 
shows that, in fact, significant similarities exist between the two plots. The modeled 
output signal displays the same jumps in its peaks and troughs as those of the 
oscilloscope output signal. Moreover, the input and output waveforms of the model 
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intersect at almost the same locations as the measured waveforms. Despite these 
parallels, one particular divergence existing between the two output signals is the fact that 
V ou t of the model never takes on both positive and negative values i.e., V ou t is either 
greater than or equal to zero or less than or equal to zero, depending on the chosen values 
of the parameters. This fact will be noted in all the test cases. Further improvements in 
the model may resolve this issue. Even though the amplitude of V ou t is not correct since 
frequency coefficient has been significantly increased for the purpose of depicting the 
details of V out , overall, V ou , of the model has the same shape as the oscilloscope’s V out and 
nearly the same values when frequency coefficient is small enough. 

Test Case II 

For the second test case, Vdd was again set to 0.709V, the amplitude of V in was decreased 
to 2V, V in.offset was set to 2V, the frequency was 1MHz, and the load resistance was 
selected to be 600 1<Q. The following table lists the values of the user-defined parameters. 

Table 3 Initial parameter values for Test Case II 


Parameter 

Value 

Vdd 

0.709V 

Rload 

600kf2 

Frequency 

1MHz 

V in, amp 

2V 

V in,offset 

2V 

Phase 

-0.75 

iDmax 

6.00E-04A 

V p,pos 

2.5V 

Vp.neg 

-0.5V 

k 

1 

Frequency Coefficient 

0.0002 


In this test case, a negative value for phase was chosen since the output signal 
lags the input signal. The magnitude of phase is much smaller for this test case than for 
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the previous test case since, as can be seen in the oscilloscope output of Figure 5, the 
input and output signals are out of phase by a small amount. The removal of the clipping 
effect was again noted when the value of k was decreased. However, decreasing k below 
1 resulted in more rounded peaks and troughs than desired, so k was set to 1 . The value 
of 0.0002 was selected for frequency coefficient because it led to the smallest V ou t,amp 
while maintaining the desired shape of V out . Figure 6 shows the modeled output. 



Figure 5 Oscilloscope output for Test Case II 


FeFET Amplifier Model 
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Figure 6 Modeled output for Test Case II based on initial parameter values 


The modeled output displays the same phase shift present in the oscilloscope 


output. Moreover, V ou t,a mp of the model is as close as possible to V ouUamp of the 
oscilloscope output, which is 74mV, while displaying the same trends as the empirical 
output signal. To further verify the accuracy of the model, plots were obtained for this 
test case when frequency coefficient was increased to 0.002. Figure 7 shows the modeled 
output, and the table below summarizes the parameter values. 

Table 4 New parameter values for Test Case II 


Parameter 

Value 

Vdd 

0.709V 

Rload 

600kf2 

Frequency 

1MHz 

V in, amp 

2V 

V in, offset 

2V 

Phase 

-0.75 

iDmax 

6.00E-04A 

V p,pos 

2.5V 

Vp, ne g 

-0.5V 

k 

1 

Frequency Coefficient 

0.002 
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Figure 7 Modeled output for Test Case II based on new parameter values 

The plots of Figure 7 look very similar to those generated by the oscilloscope. 
The modeled plots display the same shape and intersect at almost the same points as the 
oscilloscope plots. It can be concluded that the plots of Figure 6 have the correct, 
expected shape but were too small to be effectively compared to the measured results. 

Test Case III 

For the final test case, Vdd remained 0.709V, the amplitude of V in was further decreased 
to IV, 3 V of offset were applied to V in , the input frequency was 100kHz, and R/ oa d was 
560kfl These values and the other user-defined parameters are summarized in Table 5. 
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Table 5 Initial parameter values for Test Case III 


Parameter 

Value 

Vdd 

0.709V 

Rload 

560kf2 

Frequency 

100kHz 

V in, amp 

IV 

V in,offset 

3V 

Phase 

1.5 

iDmax 

6.00E-04A 

V p,pos 

2.5V 

Vp.neg 

-0.5V 

k 

2 

Frequency Coefficient 

-0.00015 


This list of parameter values is different from those presented previously in 
several aspects. Her e, frequency coefficient is a negative value. It was determined that a 
negative frequency coefficient most accurately represents the fact that the maximum 
value of V out on the oscilloscope is 75mV, and its minimum value is -50mV, as seen in 
Figure 8. In this figure, the output signal is plotted on a coordinate system with 50mV 
per division. Another point of interest is the fact that the measured output signal is 
lagging the input signal, but a positive value of phase was selected. This can be 
attributed to the fact that the negative frequency coefficient flipped the plot of V out across 
the x-axis, thereby making V out lead V in . This accounts for the positive value of phase, 
since a leading signal is given by positive phase. The parameter k was set to 2 because 
larger values for this parameter resulted in the clipping of the output signal. Thus, 
comparing the modeled plots of Figure 9 with the measured plots of Figure 8, it can be 
seen that the modeled V ou , is close in value to the measured V out . 
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Figure 8 Oscilloscope output for Test Case III 
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Figure 9 Modeled output for Test Case III based on initial parameter values 
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To obtain plots of V ou t and V in that closely resemble the oscilloscope plots, several 
changes were made to the parameter values. First, frequency coefficient was changed to a 
positive value, and its magnitude was increased to 0.0035. Since frequency coefficient 
was then positive, phase was changed to a negative value so that V out lagged V in as shown 
in the oscilloscope plots. A value of 1.3 was chosen for k because higher values resulted 
in more rounded peaks and troughs for V out than desired. The following table lists the 
new parameter values. The modeled plots are shown in Figure 10. 

Table 6 New parameter values for Test Case III 


Parameter 

Value 

Vdd 

0.709V 

Rload 

560kf2 

Frequency 

100kHz 

V in, amp 

IV 

V in, offset 

3V 

Phase 

-1.5 

iDmax 

6.00E-04A 

V p,pos 

2.5V 

Vp, ne g 

-0.5V 

k 

1.3 

Frequency Coefficient 

0.0035 
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Figure 10 Modeled output for Test Case III based on new parameter values 

These plots for V in and V out closely resemble the empirically-derived plots. As 
with the previous test cases, these plots verify the correctness of the model. It is clear 
that the desired results are obtainable using this model. 

CONCLUSION OF DATA ANALYSIS 

From these test cases, it can be concluded that the model’s results are not only 
comparable but similar to the oscilloscope outputs. The model’s output plots have the 
same shape and exhibit the same trends as those produced by the oscilloscope. 
Examining the mathematical equations behind the model brings to light their simplicity 
and the model’s ease of use. Therefore, this model is accurate enough to be used now 
and efficient enough to be easily improved for further study. 
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